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Pion form factor with twisted mass QCD 

Abdou M. Abdel-Rchim^ and Randy Lewis^ 

'^Department of Physics, University of Regina, Regina, SK, S4S 0A2, Canada 

The pion form factor is calculated using quenched twisted mass QCD with (3 — 6.0 and maximal twisting angle 
a; = ±-|. Two pion masses and several values of momentum transfer are considered. The momentum averaging 
procedure of Frezzotti and Rossi is used to reduce lattice spacing errors, and numerical results are consistent with 
the expected 0{a) improvement. 
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1. INTRODUCTION 

Lattice QCD with a twisted mass term 
(tmQCD) is a computationally efficient method 
for eliminating the exceptional configurations 
that plague light Wilson quarks, f In addi- 
tion, at maximal twist, many quantities are 
0(a) improved and others become 0{a) improved 
through a momentum averaging procedure. 
The pion form factor provides an opportunity to 
explore this improvement procedure. 

With only two valence fermions and no discon- 
nected lattice diagrams^, the pion form factor is 
an appealing laboratory for studies of the tran- 
sition between the perturbative and nonpertur- 
bative regimes of QCD. Experimental measure- 
ments are available for comparison, and current 
experiments at Jefferson Lab are exploring higher 
momentum transfers.^ 

Recent lattice studies of the pion form factor 
have considered various lattice actions, with and 
without 0{a) improvement, and the effects of 
0{a) terms are found to be non- negligible. |S] 

In this work, we present results from pion form 
factor calculations using tmQCD and compare to 
existing lattice studies. 



2. METHOD 

Computations are performed with the f3 = 6.0 
Wilson gauge action and the twisted mass ac- 
tion for a degenerate doublet of up and down 
quarks with no Symanzik improvement (clover) 
term.pP A twisting angle of ±^ is obtained by set- 



ting the hopping parameter to its critical value, 
Kc = 0.156911. PI The quark mass is then de- 
termined by the remaining parameter /i in the 
tmQCD action, and results are reported here for 
two options, 1^1 = 0.030 and 0.015, correspond- 
ing to pion masses near 660 and 470 MeV re- 
spectively. We use 100 configurations that are 
16^ X 48 with periodic boundary conditions. The 
GMRES-DR matrix inverter^, which deflates 
the smallest eigenvalues and systematically im- 
proves them upon successive restarts of the stan- 
dard GMRES iteration, was used throughout this 
work. 

The pion form factor F{Q^) is defined by 
(^+(P/)|jM(0)k+(p^O) = F{Q'){p,+pf), (1) 

where jp,{0) is a conserved vector current evalu- 
ated at the spacetime origin, pi and pf are the 
initial and final pion (Euclidean) 4-momenta re- 
spectively. Pi and Pf are the corresponding 3- 
momenta, and = (pf—pi)^ is the 4- momentum 
transfer. To compute the above matrix element 
on a spacetime lattice, one can use the three point 
correlator displayed in Fig. ^ A source with pion 
quantum numbers is placed at Xi, a sink at x/, 
and a vector current is inserted at x. The pion 
form factor is extracted from a simultaneous sin- 
gle exponential fit to the long time range of the 
two and three point correlators given by: 

G^.{U,t,p) - ^e-^(^~-^^')-P(0|^(a;)^t(^^)|o) 



cosh 



t-U--]E 



(2) 



1 



2 



quark 




anti-quark 

Figure 1 . Three point correlator for the pion form 
factor. 
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where 

(0|(/.(x)|7r+(p)) = Ze*f^ 



('r+(p/)|j^(0)k+(p.))(3) 



(4) 



4>(x) is a local interpolating field operator with 
7r+ quantum numbers, T is the temporal extent 
of the lattice and E, E^, Ef denote pion energies. 
In practice, we compute one propagator from Xi 
to X, and a double propagator from Xi to Xf to 
X. Three different options for the double prop- 
agator were studied, corresponding to the local 
pseudoscalar operator at x/ having momentum 
Pf = (0,0,0), (0, 0,pmin) and (0,0, -pmin), where 
Pmin = X (L=16). In all cases, Xi and xj are 
separated by 15 time slices. The conserved vector 
current was used at x and a local pseudoscalar at 
Xi. Smeared operators have been used routinely 
in pion form factor studies but we chose local 
operators for this first consideration of tmQCD. 

3. RESULTS 

Analyzing the pseudoscalar-pseudoscalar and 
vector- vector two point correlators, we found that 



Table 1 

Pseudoscalar and vector meson masses. 
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• Ground state energy, n=0.015 

- ■ Continuum dispersion relation, n=0.015 

■ Lattice dispersion relation 
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Figure 2. Pion dispersion relation at — 0.015. 
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Figure 3. Form factor at — 0.030 compared to 
vector meson dominance (VMD). 
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Figure 4. Form factor at |^| = 0.015 compared to pigure 5. Comparing quenched twisted mass re- 

vector meson dominance (VMD). suits with earher quenched calculations. 



a reasonable signal could be obtained for mo- 
menta < 4p^;jj. FigureElshows the dispersion 
relation for the ground state pseudoscalar meson 
at = 0.015. Table n lists results for the pseu- 
doscalar and vector meson masses in lattice units. 

To extract the form factor at higher Q^, 
nonzero sink momentum is important, but then 
improvement requires momentum averaging over 
positive and negative momenta. f^I Figures |31 and 
^show our results for the pion form factor as a 
function of Q^. 

A comparison with existing literature is shown 
in Fig. Notice that unimproved Wilson re- 
sults are systematically below experiment (which 
follows vector meson dominance in this region of 
Q^)- Decreasing decreases F{Q^) even fur- 
ther. However, the tmQCD results are consistent 
with experiment and with 0(a) improved actions. 
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